Assessment of a quantum phase gate operation based on nonlinear optics 
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We analyze in detail the proposal for a two-qubit gate for travelling single-photon qubits recently 
presented by C. Ottaviani et al. [Phys. Rev. A 73, 010301(R) (2006)]. The scheme is based on an 
ensemble of five- level atoms coupled to two quantum and two classical light fields. The two quantum 
fields undergo cross-phase modulation induced by electromagnetically induced transparency. The 
performance of this two-qubit quantum phase gate for travelling single-photon qubits is thoroughly 
examined in the steady-state and transient regimes, by means of a full quantum treatment of the 
system dynamics. In the steady-state regime, we find a general trade-off between the size of the 
conditional phase shift and the fidelity of the gate operation. However, this trade-off can be bypassed 
in the transient regime, where a satisfactory gate operation is found to be possible, significantly 
reducing the gate operation time. 

PACS numbers: 03.67.Mn, 42.50Gy, 42.65.-k 



I. INTRODUCTION 

A promising system for quantum computation is to use 
single photons to encode quantum information [JJ. This 
is due to the photon's robustness against decoherence and 
the availability of single-qubit operations. However, it 
is difficult to realize the necessary two-qubit operations 
since the physical interaction between photons is very 
small. Linear optics quantum computation and non- 
linear optical processes involving few photons have been 
proposed to circumvent this problem. The first is a prob- 
abilistic scheme implicitly based on the nonlinearity hid- 
den in single-photon detectors, while the second is based 
on the enhancement of photon-photon interaction either 
in cavity QED configurations ,3| or in dense atomic me- 
dia exhibiting electromagnetically induced transparency 
(EIT) 4J. In this latter case, optical nonlinearities can 
be produced when EIT is disturbed, either by introduc- 
ing additional energy level(s) H0, or by mismatching 
the probe and control field frequencies [all . 

The scope of this paper is to assess the performance 
of a two-qubit quantum phase gate (QPG) for travelling 
single photon qubits [l(J, HH D3 UH , based on the 
cross-Kerr nonlinearity which is generated in a five-level 
atomic medium. In a QPG, one qubit gets a phase con- 
ditional to the other qubit state according to the trans- 
formation |«) 1 1 j) 2 — > exp{«0jj} |i)i|j) 2 where 
= 0,1 denote the logical qubit bases. This gate 
is universal when the conditional phase shift (CPS) 



(1) 



is nonzero, and it is equivalent to a CNOT gat e up 
to local unitary transformations when <j) = it [la, Hq |. 
Most of the literature focused only on the evaluation 
of the CPS and on the best conditions for achieving 
= 7T E3, E3, E3, E3, E3 , while the gate fidelity, which is 
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the main quantity for estimating the efficiency of a gate, 
has been evaluated in the full quantum limit in Ref. 0] 
for the first time. Here we provide the details of the cal- 
culation of the fidelity and the CPS of Ref. [H, which 
showed the presence of a general trade-off between a 
large CPS and a gate fidelity close to one, hindering the 
QPG operation, in the stationary state. However, we 
shall see that this trade-off can be partially bypassed in 
the transient regime, which has never been considered 
before in EIT situations, still allowing a satisfactory gate 
performance. 

The qubits are given by polarized single-photon wave 
packets with different frequencies, and the phase shifts 
(j>ij are generated when these two pulses cross an atomic 
ensemble in a five- level M configuration (see Fig.^). The 
population is assumed to be initially in the ground state 
|3). From this ground state, it could be excited by either 
the single-photon probe field, with central frequency lu p 
and coupling to transition |3) <-> |2), or by the single- 
photon trigger field, with central frequency u) t and cou- 
pling to transition |3) <-> |4). We assume that the five 
levels are Zeeman sub-levels of an alkali atom, and that 
both pulses have a sufficiently narrow bandwidth. In 
this way, the Zeeman splittings can be chosen so that the 
atomic medium is coupled only to a given circular po- 
larization of either the probe or trigger field, while it is 
transparent for the orthogonally polarized mode, which 
crosses the gas undisturbed [Tj. As a consequence, the 
logical basis for each qubit practically coincides with the 
two lowest Fock states of the mode with the "right" po- 
larization, and |lj) (j — p,t), while the "wrong" 
polarization modes will be neglected from now on. 

A classical pump field, with frequency oj\ and Rabi 
frequency f2i, couples to the transition |1) <-> |2), while 
a second classical pump field, with frequency lo^ and 
Rabi frequency SI4, couples to the transition |4) |5) 
(see Fig.^). We consider a cylindrical, quasi-lD, atomic 
medium with the two classical pump beams propagat- 
ing along its axis, collinear with the two quantum fields 
in order to avoid Doppler broadening. When the probe 
field is on two-photon resonance with the pump field with 
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FIG. 1: Energy levels of the M-scheme. flj are the Rabi fre- 
quencies of classical fields, while g Pt t denote couplings of the 
quantized probe and trigger fields to their respective transi- 
tions. Sj are the detuning of the fields from resonance. 

Rabi frequency Hi, and the trigger field is on two-photon 
resonance with the pump field with Rabi frequency Q4, 
the system exhibits EIT for probe and trigger simultane- 
ously. This simultaneous EIT condition is achieved when 



Si = 82, S3 = S4, (2) 

where the detunings Sj are defined by 

E 2 -Ex = hu>i (3a) 

E 2 -E 3 = huj p + hS 2 , (3b) 

£4 -£3 = huj t + nS 3 , (3c) 

E4 — E§ = T1UJ4 + TlSa. (3d) 



A nonzero CPS occurs whenever a nonlinear cross-phase 
modulation (XPM) between probe and trigger is present. 
This cross-Kerr interaction takes place if the two-photon 
resonance condition is violated. For small frequency mis- 
match 612 = S\ — 82 and £34 = S3 — 84 (both chosen to be 
within the EIT window), absorption remains negligible 
and the cross-Kerr interaction between probe and trig- 
ger photons may be strong. The consequent CPS may 
become large, of the order of tt, if the probe and trigger 
pulse interact for a sufficiently long time. If the two single 
photon pulses enter simultaneously the atomic medium, 
their interaction time U nt is optimized when the group 
velocities of the two pulses are equal, so that tint — L/v g , 
where v g is the common group velocity of the pulses and 
L is the length of the gas cell. The inherent symmetry 
of the scheme guarantees perfect group velocity match- 
ing for probe and trigger whenever 8\ = 84, 82 = 83 and 
g p /fli = gt/Qi, where gj = /j,j^/ojj/2he Vj (j = p,t) is 
the coupling constant between the quantum mode with 
frequency ujj and mode volume Vj , and the corresponding 
transition with electric dipole moment fj,j. 

The importance of group velocity matching for achiev- 
ing a significant nonlinear phase shift was first pointed 
out in Ref. Ifj], which suggested to use a mixture of two 
different atomic species to achieve this goal. The first 
kind of atoms generates XPM by means of a four-level N 
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FIG. 2: (Color online) A schematic plot of the assumed single- 
photon pulse propagation through the gas cell of length L and 
diameter d. The pulse length is assumed to coincide with the 
cell length L, and the pulse waist w is assumed to be of order 
of cell diameter d. 

scheme '0\ , in which however only the probe field under- 
goes EIT and it is slowed down, while the second kind of 
atoms realizes a three-level A scheme able to slow down 
the trigger pulse. Group velocity matching is achieved by 
means of an accurate but difficult control of the atomic 
densities. A different way of achieving group velocity 
matching, but which is still asymmetric for probe and 
trigger has been proposed in Ref. [TTj . which employs a 
five-level M scheme similar to the one discussed here, but 
with a different atomic population distribution. In that 
case, the two group velocities can be tuned and made 
equal simply by tuning the frequencies and intensities of 
the two classical pump fields. Instead, Ref. consid- 
ered a six level scheme in which probe and trigger are 
affected by EIT and XPM in a symmetric fashion, so 
that the corresponding group velocities are equal by con- 
struction. The present proposal achieves group velocity 
matching just in the same way (see also Refs. where 
a four-level tripod system, symmetric between probe and 
trigger, has been proposed for XPM). 

The paper is organized as follows. In Sec. [n] we de- 
scribe the model used in the remainder of the paper. 
Sec. 11111 shows the results of a perturbative calculation 
for the CPS. These are used as a motivation to pass to a 
density matrix based calculations in Sec. IIVI describing 
a QPG operation in a steady-state. Then, the transient 
regime in explored in Sec. [V] while in Sec. IVII a scheme 
for the experimental verification of the QPG operation is 
discussed in detail. Conclusions are given in Sec. IVIII 

II. MODEL 

In this Section, we present the model we have adopted 
for a full quantum description of the interaction of the 
two single-photon wave-packets with the atomic medium 
possessing the level structure outlined in Fig. ^ To this 
end, we make the following two assumptions which, even 
though not simple to realize experimentally, are more 
technical than physical in nature: 

1. We assume perfect spatial mode matching between 
the input single-photon pulses entering the gas cell 
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and the optical modes naturally excited by the 
driven atomic medium, and which arc determined 
by the geometrical configuration of the gas cell and 
of the pump beams • This allows us to describe 
the probe and trigger fields with the right polar- 
ization in terms of single travelling optical modes, 
with annihilation operators a Pt t- 

2. We assume that the pulses are tailored in such a 
way that they simultaneously enter gas cell and 
completely overlap with it during the interaction 
(see Fig. [5J. This means that their length (com- 
pressed inside a medium due to group velocity re- 
duction) is of the order of the cell length L and 
their beam waist is of the order of the cell radius. 
In this way, the two pulses interact with all N a 
atoms in the cell at once, and moreover, one can 
ignore spatial aspects of pulse propagation. 

With these assumptions, and neglecting dipole-dipole in- 
teractions, the interaction picture Hamiltonian may be 
written as 

H = hei 2 Sn + hS 2 S 22 + M3S44 + he 3i S 55 (4) 

+nn 1 t/K (£21 + £12) + frg P \/K (a p s 2 3 + s 32 a P ) 

+hgtVNa (a t S ia + S 3i afj + hQ iV ^Na [Si5 + £54) , 
where we have defined the collective atomic operators 
1 N ° 

Sh = —=^2<J l kU k^l = l,...,5, (5a) 
V iVa i=l 
N a 

Skk = J2 a kk' ( 5b ) 

i 

with ah = being the operator switching between 

states k and I of the ith atom. The initial state of the sys- 
tem corresponds to a probe and a trigger single-photon 
pulse with generic polarization, simultaneously entering 
the medium in which all the atoms are initially in state 
1 3). Since we consider only the polarization mode in- 
teracting with the medium, both for the probe and the 
trigger, the initial state can be written as 

Win) = (g> 13> ( ® (coo|O p > ® |0t> + cdi|0„> (81 ]l t > 

i=l 

+eio|lp)®|Oi)+cu|lj,)® |l t )). (6) 

Due to the above assumptions, the passage of the two 
pulses through the atomic medium of length L corre- 
sponds to the time evolution of this state, for a time 
Unt = L/vg, according to the master equation [5| 

N a 

p = -%l H > d + X] I? £ ( 2a ikP a ik - a lkp - P a lk) . 

k j=l 

+EirE ( 2 4p4I - °%4ip KM*) . (?) 

kl ] = 1 



including not only the coherent interaction described by 
the Hamiltonian of Eq. J2J, but also the spontaneous 
emission from the excited states I = 2, 4 to the ground 
states k = 1,3,5 (jki denotes the corresponding decay 
rate) and the dephasing of levels |fc), k — 1, 2, 4, 5, with 
dephasing rate jkk- Typically the dephasing rates are 
much smaller than the decay rates, jkl 3> 7fcfc, V k, I. 

Since the initial state of Eq. © contains at most two 
excitations, the coherent time evolution driven by Eq. Q 
is simple and restricted to a finite-dimensional Hilbcrt 
space involving few symmetric collective atomic states. 
In fact, each component of the initial state of Eq. © 
evolves independently in a different subspace. The com- 
ponent with no photon is an eigenstate of H and does not 
evolve. The ®^ |3) i |0 P ) £g> |l t ) component evolves in a 
three-dimensional Hilbert space which it spans together 
with the two states le^ ' ^) and (eg 0,0 ''). Here, we have 
defined the symmetric collective states 

1 N « 

\ e (n P ,n t)) = J2\^2,...,ri,...,3 Na )^\n p )^\TH), 

where r — 1,2,4,5. In a similar fashion, the compo- 
nent with only one probe photon evolves in a three- 
dimensional Hilbert space spanned by the three states 

I 3 ), lip) ® |0t), |4°' 0) ) and |4°' 0) ). The component 
with one probe and one trigger photon evolves in the 
five dimensional subspace spanned by the four collective 
states |ef '^), le^ ' 1 ^), le^ 1,0 ^), and leg 1 ' ^) and the state 

<g&|3)Jl P )®|l t )- 

Decoherence effects, and more specifically spontaneous 
emission from each atom complicates this dynamics. 
However, we are in the weak excitation limit where, for 
I 7^ 3, {(T°n) ~ N^ 1 < 1, as shown by the fact that 
the Hamiltonian dynamics involve only the symmetric 
atomic states of the form of Eq. JSJ. This limit al- 
lows a drastic simplification of the effective time evo- 
lution. Following Duan et al. JJj, we can introduce 
Fourier transforms of the individual atomic operators 

4l = E&~Ml e<jWJVa /V2^, where s li = Ski are the 
collective operators defined in Eq. (|5a|l . The sum over 
the atoms in Eq. then transforms to the sum over the 
collective atomic modes with index /1. In the weak ex- 
citation limit, the operators s kl approximately commute 
with each other. This means that they represent inde- 
pendent collective atomic modes, and one can trace over 
the fi ^ modes, so that the spontaneous emission term 
in the master equation becomes 

X -y ( 2S kipSti ~ S^Skip - pS^SkiJ , (9) 

kl 

where the sum is now over the six "collective" spon- 
taneous decay channels only, each characterized by the 
single-atom decay rate 7/.;. A similar argument applies 
to the dephasing term in the master equation Q. In fact, 
if we restrict to the subspace of the symmetric collective 
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states of Eq. (|SJ involving only single atomic excitations, 
we can approximate in the dephasing terms of the master 
equation, 



E 



ikk J2 a 



j j 

kkP a kk 



^ IkkSkkpSkk, 



(10) 



is obtained by using the fact that the phase shifts (pij of 
Eq. (JIJ are given by combinations of the phases of the 
off-diagonal matrix elements (in the Fock basis) of the 
reduced density matrix of the probe and trigger modes, 
Pf(t) = Tr atoms {p(t)}. 

The gate fidelity is given by |16| 



where Skk is given by Eq. (|5bjl . Using Eqs. © and ffTTTfl . 

the master equation of Eq. Q in the weak excitation 
limit becomes 



= -- [H, p] + Yl "IT V^SkkpSkk - SkkP - pSkk 



+ Ey (^mpsI - slSkw - pslsu) , (ii) 

kl 

that is, it involves only the operators of the collective 
atomic mode with index p = 0. This actually means 
that the single photon probe and trigger pulses excite 
only a restricted number of collective atomic states, so 
that the atomic medium behaves as an effective single 
5-level atom, with a collectively enhanced coupling with 
the optical modes gj \J N a , but with single- atom decay 
rates "fki, dephasing rates "fkk, Rabi frequencies Qi, and 
detunings Si. 

Spontaneous emission causes the four independent 
Hilbert subspaces corresponding to the four initial state 
components to become coupled. Moreover, the "cross" 
decay channels |4) — > |1) and |2) — > |5) couple the above- 
mentioned collective states with six new states, le^ 1 ' ^), 
(populated if 741 7^ 0), and leg ' 1 - 1 ), 

e£ Q,1) ), |4°' 2) ) (populated if 725 ^ 0). Therefore Eq. JTTJ 
actually describes dynamics in a Hilbert space of dimen- 
sion 18, which we have numerically solved in order to 
establish the performance of the QPG. Notice that, due 
to the combined action of the cross-decay channels and 
of the Hamiltonian QJ, the states |e^' ^), |e 2 ), leg ' 1 ^) 

and le^ 0,1 "*) are coupled also to doubly excited atomic col- 
lective states without photons which are neglected by our 
treatment. However, as we shall see below in the paper, 
a good QPG performance is possible only when spon- 
taneous emission events are rare. Under this condition, 
the probability to populate these doubly excited atomic 
collective states during the atom-field interaction is com- 
pletely negligible, and therefore our model based on the 
effective single five-level atom description provided by 
Eq. i|ll|) is essentially correct. 

In summary, the model outlined above relies on the 
single-photon nature of the excitations. In this case, 
the collective operators (0 effectively switch between 
the states making the superposition \ipi n ) ©, and the 



4°' 2) > 



symmetric collective states 



(n p ,n t ) 



©. This is central 



to the reasoning leading to the effective master equa- 
tion itTT|) . 

To characterize the QPG operation, we calculate the 
CPS (j) °f Eq. and the fidelity of the gate. The accu- 
mulated CPS 4> as a function of the interaction time tint 



where 



^(*) = V^«(*)IM*)I^«(*)>, (12) 



\i>id{t)) = c O oexp{i0oo(imt)}|Op,O t ) 
+c i exp{i(/>oi(£)}|0 p , l t ) 
+C10 exp{^ 10 (i)}|l p ,0 t ) 
+c 11 exp{i(/) ll (t)}\l p ,l t ) (13) 

is the ideally evolved state from the initial condition 
© , with phases faj (t) evaluated from p f (t) as discussed 
above. The overbar denotes the average over all initial 
states (i.e., over the cy, see Poyatos et al. 20}). The 
above fidelity characterizes the performance of the QPG 
as a deterministic gate. However, one could also consider 
the QPG as a probabilistic gate, whose operation is con- 
sidered only when the number of output photons is equal 
to the number of input photons. The performance of 
this probabilistic QPG could be experimentally studied 
by performing a conditional detection of the phase shifts, 
and it is characterized by the conditional fidelity !F c (t), 
which will be discussed in Sec. IIVI 



III. PERTURB ATIVE REGIME 

The conditional fidelity is always larger than the un- 
conditional one, but they become equal (and both ap- 
proach 1) for an ideal QPG in which the number of pho- 
tons is conserved and all the atoms remain in state |3). 
This ideal condition is verified in the limit of large de- 
tunings Sj 3> 7fej so that spontaneous emission is signif- 
icantly suppressed and can be neglected, and very small 
couplings gj V N a <C flj- In this limit, each component 
of the initial state of Eq. JfJJ practically coincides with 
the dark state of the four independent Hamiltonian dy- 
namics discussed in Sec. [H] The system with the initial 
state containing zero probe and trigger photons does not 
evolve, i.e. stays in the initial state i= ?j |3)^ |0 P ) (g> |0 t ) 
all the time. The subsystem containing one probe and 
zero trigger photons as the initial state evolves accord- 
ing to a reduced three-dimensional Hamiltonian, which 
in the basis formed by the states 0^ |3) i |l p ) <g> |0 t ), 



,(0,0) 



) and 



3 (0,0)v 



is given by 



H p = I g p ^/Na 82 



0_ g P VK 

tti £12 



(14a) 



Similarly, the subsystem containing one trigger and zero 
probe photons as the initial state evolves according to 



K 




FIG. 3: Conditional phase shift (left) and the fidelity of Eq. 1121 (right) as a function of the interaction time for N a = 10 B , 
Si = S3 — 7.57, ei2 = £34 = 0.057, 9p — 9t = O.OOII7, Q,i = 0,4 = 1.8757 an d 7fcfc = Jph = 10~ 3 7, V k. We have taken equal 
decay rates, 721 = 723 = 725 = 741 = 743 = 745 = 7/3, with 7 = 2n x 6 MHz. Left: solid line represents the phase shift, as 
calculated from the full master equation, while the dashed line gives the perturbative prediction of Eq. 1161 . Right: solid line 
is the unconditional gate fidelity J-(t), while the dot-dashed line is the conditional one, J~ c (t). 



a reduced three-dimensional Hamiltonian, which in the 



basis formed by the states i= ?i 13)^0 



11/ 



and 



,(0,0) v 



is given by 




,(o,o) v 



(14b) 



Finally, the subsystem containing initially one photon 
each in probe and trigger modes evolves according to a 
reduced five-dimensional Hamiltonian, which in the basis 
formed by the states \e( 



(0,1) 



,(°,1)\ 



,(1,0)* 



, and 



(1,0), . . , 

25 )) is s iven b y 



(SCiISUp/ 



1, 
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£34 



(14c) 



The phase accumulation experienced by the various com- 
ponents of the quantum state of the fields will be pro- 
portional to the eigenvalues of these matrices. The four 
phase shifts 4>ij can be evaluated as a fourth-order per- 
turbation expansion of the eigenvalue corresponding to 
the dark state in each subspace, multiplied by the inter- 
action time tint- The CPS is then calculated as 

<t> = ( x H pt - A ffp - X Ht )Unt, (15) 

where the A's denote the eigenvalues of the correspond- 
ing reduced Hamiltonian, with \H pt 0n, ^h p *-> <fiw, 
A.ff t <-> <poi and 0oo = 0, m agreement the with general 
definition of Eq. . Following this procedure results in 
the following CPS 



gWtNlt 



ini 



(^34^3 



e 34(£i2 



fll)(e 12 6 2 

-n?) , 



-n?) 

£12 (e 2 



34 



Ql) 



(£12^ - n?) (e 3 4^3 - n?) 



(16) 



This prediction is verified by the numerical solution of 
Eq. (|llfl in the limit of large detunings and small cou- 
plings. However the resulting CPS is too small, even for 
very long interaction times (i.e., long gas cells): for exam- 
ple, for g p ,tVNa = 0.5 MHz, ei 2 , 34 = 1.9 MHz, fi li4 = 65 
MHz and 62. 3 =1.9 GHz, we obtain a tiny CPS of only 
3 x 1CP 4 radians when tint = 10~ 4 s, which corresponds 
to L ~ 30 km. This is not surprising because this limit 
corresponds to a dispersive regime far from EIT. In this 
regime, transparency is achieved by means of a strong 
coupling field, producing a well-separated Autler-Towncs 
doublet 5j. At the same time, the size of nonlinearity is 
small due to the extremely weak coupling of the quan- 
tized fields to their respective transitions. The results 
are, therefore, not different from those expected from 
XPM in a standard optical fiber |2l]]. Therefore, one 
has to explore the non-perturbative regime of larger cou- 
plings in order to exploit the low-noise, large-nonlinearity 
properties of EIT and achieve a satisfactory QPG oper- 
ation. 



IV. STEADY-STATE QPG OPERATION 

A large amount of work exploring EIT-based nonlinear 
optical phenomena considers the steady-state of a generic 
EIT-based system as being the natural state in which to 
predict and test different phenomena 0, S IS EH • We 
shall see in this Section that it is not possible to achieve 
a satisfactory QPG performance in such a steady-state 
regime. 

In this Section, we analyze the performance of the 
QPG at the steady state. To this end, we consider two 
different parameter regimes: (i) the regime of long in- 
teraction times, a natural extension of the perturbation 
analysis of Sec. IIIII and (ii) the regime of short interac- 
tion times, corresponding to a non-perturbative regime 
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FIG. 4: Conditional phase shift (left) and the fidelity (right) as a function of the interaction time for N a — 10 s , Si = S3 = 9.57, 
£12 = £34 = O.27, g p — g t — O.OI87 and J7i = O4 = 197 and ^ kk — j p h = 10 _3 7, V k. We have taken equal decay rates, 
721 = 723 = 725 = 741 = 743 = 745 = 7/3, with 7 = 2ir x 6 MHz. Left: solid line represents the phase shift, as calculated from 
a density matrix, while the dashed line gives the eigenvalue solution. Right: solid line is the unconditional gate fidelity T(t), 
while the dot-dashed line is the conditional one, J~ c (t). 



with strong atom-field coupling. 

A. Long Interaction Time 

Naturally extending the perturbative analysis, we solve 
the master equation (|llf> . and show the results in Fig. [3] 
Fig. |3 (left) shows the result for the conditional phase 
shift. Solid line has been calculated from the solution 
of Eq. dJ, as explained in Sec. |HJ The dashed line is 
the 'benchmark' solution, obtained from the eigenvalues 
of the associated Hamiltonian, by using Eq. I|15p. The 
eigenvalues of the Hamiltonians of Eqs. (|14|) have been 
calculated numerically for the set of parameters of Fig. [3] 

It is evident that the 'benchmark' solutions offers a rea- 
sonably good estimate for the size of the CPS. The exact 
dynamics driven by the master equation l|llfl presents 
an additional oscillatory behavior both on a short time 
scale (transient processes), and on a long-time scale. The 
longer time scale comes from the fact that to induce the 
cross-Kerr nonlinearity, one has to detune the fields away 
from the dark resonance. This detuning is very small and 
is seen in the oscillations on a long time-scale. As both 
probe and trigger fields are detuned by the same amount, 
only one frequency of long-time oscillations is observed. 

In Fig. |3 (right), fidelities (averaged over all possible 
two-qubit initial states) are shown in two cases. Both are 
calculated by using Eq. Q12|). but they differ in the way 
Pf(t) is defined. The solid line in Fig. [3] (right) refers to 
the unconditional fidelity F(t), which is calculated from 
Eq. {nj by taking p f (t) = Tr at {p(t)}, where p(t) is the 
solution of the master equation l|llf> . The unconditional 
fidelity quantifies the performance of the QPG as a de- 
terministic gate for single-photon qubits. 

The dot-dashed line in Fig. [3] (right) refers to the 
conditional fidelity T c (t), which is evaluated accord- 
ing to Eq. (|T2)l . but with Pf(t) replaced by py(t) — 

^at{\^nj(t)){^nj(t)\}/(^nj{t)\^nj(t)), where \4> n j(t)) IS 



the (non-normalized) evolved atom-field state condi- 
tioned to the detection of no quantum jumps [23|, i.e., 
of no photon loss by spontaneous emission. This fidelity 
can be measured by post-selecting those measurement re- 
sults that conserve photon number, i.e., discarding those 
data sets where at least a photon from the initial two- 
qubit state has been lost to the environment. The condi- 
tional fidelity quantifies the performance of the QPG as 
a probabilistic two-qubit gate. 

In Fig. 13 we have found at best a CPS of ~ it in 
correspondence with fidelities T(U n t) and F c (U n t) equal 
to 60% and 80%, respectively. This is due to the gen- 
eral presence of a trade-off between the size of the CPS 
and of the gate fidelity ; as well as to the atomic dephas- 
ing |28| . This is an important result of our paper, which 
actually holds true in any EIT-based nonlinear optics 
systems. In fact, both the conditional and the uncondi- 
tional gate fidelity approach 1 in the limit of very small 
gj, but this limit yields a CPS which becomes appreciable 
only for unrealistically long gas cells. For the parameters 
of Fig |3 v g — 5.8 x 10 6 , which requires a cell of length 
Vg x t int = 5.8 x 10 6 m/s x (5OO/7) ~ 81 m. Therefore a 
larger CPS requires a larger ratio gjy/N a /Flj. This con- 
dition however increases the population of the collective 

atomic states |e$" p '™ t ' ) ) and leg'"*') at the expense of the 
initial atomic state |3), thus unavoidably decreasing the 
gate fidelity. Similar conclusions hold for other options, 
such as increased detunings 5j, or adjusting two-photon 
detunings e^. Therefore, just the pure coherent unitary 
evolution of the system, governed by the Hamiltonian of 
Eq. J2J causes this inherent trade-off. 

B. Short Interaction Time 

To further illustrate our findings, we calculate the 
CPS and the gate fidelities in the range of parameters 
where the total interaction time is an order of magni- 
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FIG. 5: Conditional phase shift (left) and the fidelity (right) of the QPG operation for N a — 10 6 , Si = S3 — 15j, erz = £34 = 
O.OI7, g p = gt — 0.00227, Oi = Q.4 = 47 and jkk = J v h = 10~ 3 7, V k. We have taken equal decay rates, 721 = 723 = 725 = 
741 = 743 = 745 = 7/3, with 7 = 2-7T x 6 MHz. Right: the unconditional fidelity (solid) and conditional fidelity (dashed) are 
shown. See text for details. 



tude smaller than in Sec. IIV Al The CPS and the gate 
fidelities are calculated as described in Sec. IIV Al and the 
results are shown in Fig. 0] To obtain a CPS of the or- 
der of 7T in a shorter interaction time (tint ~ 5O/7), we 
have assumed a larger ratio gj\/N a /£lj. The trade-off be- 
tween the amount of accumulated nonlinear phase shift 
and the gate fidelity is now even more pronounced: we 
find at best a CPS of ~ tt in correspondence with fidelities 
f(tint) and T c (ti n t) equal to 65% and 73%, respectively. 
As expected, having a stronger atom-field coupling en- 
hances the processes lowering the fidelity. The system 
ends up with a large CPS faster, but this is achieved with 
a final state in which the probability of loosing the probe 
and trigger photons by spontaneous emission or within 
the atomic medium is no more negligible. We notice that 
in this short interaction time case dephasing does not 
have an appreciable effect, that is, the results without 
dephasing are indistinguishable from those with dephas- 
ing shown in Fig. ^ due to fact that dephasing rates in 
a dilute gas are typically much smaller than decay rates. 
The only possible way to circumvent this trade-off is to 
explore the transient regime, which will be discussed in 
the following Section. 

V. QPG OPERATION IN TRANSIENT REGIME 

In Sec. IIVI we have found that the QPG operation of 
EIT-based nonlinear system in a steady-state is plagued 
by the trade-off between the phase shift size and the gate 
fidelity. In an attempt to find favorable conditions for the 
QPG operation, we consider the transient regime, when 
jtint < 1. As discussed above, in order to accumulate a 
significant CPS in such a short time one has to consider 
the strong coupling regime with a large ratio g^sj 'N a /flj. 
Therefore, the trade-off between fidelity and a large non- 
linear interaction is present also in the transient regime. 
However, when gjy/N a /Clj is large, the transient dynam- 
ics is characterized by Rabi-like oscillations of the atomic 



populations and of the photon number, determining, as 
a consequence, coherent oscillations of the gate fidelity. 
In such a case one cannot exclude the existence of spe- 
cial values of the interaction time U n t corresponding to 
a maximum of the gate fidelity close to one, and at the 
same time, to a value of the CPS close to ir. 

We show that this fact is actually possible in Fig. 
where we see that a CPS of ~ 7r radians is obtained in the 
transient regime for i;„ t ss O.4/7 ~ 10 ns. At the same 
interaction time, the unconditional gate fidelity (Fig. [5] 
right, full line) is about 94%, while the conditional gate 
fidelity reaches the value of 99% (Fig. [S] right, dashed 
line). The conditional gate fidelity is obtained in corre- 
spondence with a success probability of the gate equal to 
0.94, calculated from the norm of the Monte-Carlo wave 
function |23| . The probe and trigger group velocities are 
calculated to be v g ~ 3 x 10 6 ms _1 , yielding a gas cell 
length L = Vgti n t ~ 3.1 cm. The value of gj yields an 
interaction volume V ~ 2 • 10~ 3 cm 3 , corresponding to 
a gas cell diameter of about 330 /Ltm and to an atomic 
density N a /V ~ 5 • 10 10 cm" 3 . 

To give a further example and to deepen our discus- 
sion of QPG performance and the CPS-fidelity trade-off, 
we now show the optimal results for the experimentally 
available pulses produced in the experiment of Darquie 
et al. [52] (see Fig. EJ). The length of a pulse produced 
in [23| is 26 ns (~ 7" 1 )- The optimal parameters of 
Fig. give the unconditional gate fidelity (Fig. right, 
solid line) of 85%, while the conditional fidelity reaches 
89% (Fig. |B] right, dot-dashed line). The probe and trig- 
ger group velocities are calculated to be v g ~ 1.51 x 10 6 
ms^ 1 , yielding a gas cell length L = v g ti nt ~ 3.82 cm. 
The value of gj yields an interaction volume V ~ 10 _1 
cm 3 , corresponding to a gas cell diameter of about 910 
fim and to an atomic density N a /V ~ 10 9 cm" 3 . So, tak- 
ing the longer wavepackets (i.e. longer interaction times) 
means that the CPS-fidelity trade-off becomes increas- 
ingly important, and the top value of fidelity decreases 
slightly with respect to its optimal value. 
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FIG. 6: Conditional phase shift (left) and the fidelity (right) of a QPG operation for N a — 10 , #2 = <?3 = 67, £12 = £34 = 0.057, 
9p = <?« = 0.00097 an( l fh = ^4 = 7- We have taken equal decay rates, 721 = 723 = 725 = 741 = 743 = 745 = 7/3. For 
87 Rb, 7 = 2-7T x 6 MHz, giving the interaction time (i.e. pulse length) of ~ 25 ns for \<j>\ ~ w. Right: solid line denotes the 
deterministic fidelity, while dot-dashed line denotes the conditional fidelity. 



Note how in both cases (Figs. |S] and |Gj high values of 
fidelity could be obtained for a values of CPS lower than 
7r radians. This implies the possibility of the implemen- 
tation of a universal quantum gate |lq |. which requires 
only 4> 7^ 0. 

A comment about this calculation of the common 
group velocity of the two wave-packets, v g , is in order. As 
mentioned earlier, EIT is stationary phenomenon, and in 
fact, the conventional v g is a steady-state quantity which 
it is obtained from the susceptibility x according to 



(17) 



(u>o is the central frequency of wave-packet), where the 
susceptibility of the j-th field, \j U = Pi *)j is evaluated 
from the associated steady-state atomic coherence pj s as 



velocity is a signature of a sort of "non-stationary" EIT 
process. 

In order to verify that the non-stationary dynamics is 
really reminiscent of EIT, in the next subsection we com- 
pare these results with a numerical study of the three- 
level ladder atomic scheme (see Fig. 0), yielding XPM 
without EIT. Here we anticipate that we have found a 
smaller gate fidelity (~ 78%) for a corresponding set of 
parameters, providing therefore further support to the 
presence of a moderate, non-stationary EIT process in 
the transient dynamics of our five-level M scheme. This 
is because in the ladder scheme, the increase in the size 
of the nonlinearity (and thus the CPS) is accompanied 
by an increase of losses, unlike the case of any EIT-based 
scheme 6]. For this reason we also expect that our M 
scheme outperforms ladder scheme in the case of fully 
optimized fidelity for 7r phase shift. 



Xl ~ vn^T^ ■ (18) 

Instead, the above results are obtained in the transient 
regime where jtint < 1, and for this reason we have 
estimated the group velocity in a different way. We 
have evaluated the relevant time-dependent atomic co- 
herence pj (t) and the corresponding "instantaneous sus- 
ceptibility" Xj [t) from the reduced atomic density matrix 
Pred{t) = Tr fields{p(t)}, with p(t) being the solution of 
Eq. i|llfl . The corresponding "instantaneous" group ve- 
locity v g (t) has been then averaged over the time inter- 
val between and tj n t, providing in this way our esti- 
mate of the "transient" non-stationary group velocity of 
the single-photon wave-packets. For the parameters of 
Fig. this non-stationary v g is approximately equal to 
c/100 and it is about one order of magnitude smaller than 
the conventional v g obtained from the steady-state sus- 
ceptibility. This appreciable slowing down of the group 



A. The Conventional Three-Level Scheme 

The atomic ladder scheme (see Fig. 0| is well-known 
to exhibit XPM of the two fields involved @ . In order to 
achieve a reasonable size of cross-phase shift, the detun- 
ing of the intermediate state 5 P needs to be large. This 
minimizes spontaneous emission (~ but also the 

size of XPM ~ l/Sp. 

In order to evaluate the XPM in a manner comparable 
to what we have done for the M-scheme, we make similar 
assumptions and arrive at a description analogous to the 
one described in Sec. [H] The Hamiltonian is now given 
by 

H 3 = hS p S 2 2 + H$p - $t)Sa3 + ?ig P \/~Na (a P S 2 i + Swa^J 
+hg t y r K (a t S 2 3 + S 32 aty . (19) 
Following the same reasoning as in Sec. ^ we arrive at 
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the effective master equation (we neglect here atomic de- 
phasing) 

i 

p = C 3 p = -- [H 3 , p] 

+ -^- ^2S'i2pS'21 — S21S12P — pS , 2lS'l2^ 

+ "7T ( 2 'S'23P<S'32 - S23S32P - p§23§32J ) (20) 

where 712 and 723 denote the spontaneous emission rates 
from levels |2)and|3) to levels |l)and|2) respectively, 
and the operators S%j denote collective atomic operators, 
in the spirit of Eq. J5aj. 

The results of the calculation of unconditional quan- 
tities are shown in Fig. |H| for a parameter regime com- 
parable to that discussed in relation to Figs.[S] (see also 
Ref. ^tJ) f° r the five- level M scheme. The atoms are as- 
sumed to be in the collective state <S)j=i |2)j initially, as 
this is found to give better results. The reason is a more 
efficient photon-photon interaction since the initial state 
is symmetric with respect to probe and trigger photons. 
At the interaction time ti nt ~ O.I2/7, the CPS reaches 
the value ~ ir and at the same time the unconditional fi- 
delity is found to be ~ 78%. It is possible to see that the 
conditional fidelity, even though higher, always remains 
significantly lower than that obtained in the M-scheme. 

Therefore, we conclude that the optimal results for 
the QPG operation can be found in the "transient EIT- 
regime" . The general trade-off between the nonlinear 
phase shift and the fidelity is still present, but it is com- 
pensated by the transient oscillations in the populations 
of atomic levels. In fact, the numerical results show that, 
in the parameter regime under consideration, the popula- 
tion of the excited states \e 2 np '"'•') and \e£" p,nt ^) is always 
negligible, and one has coherent oscillations of the pop- 
ulation between the states \e{ v,nt ^) 1 0^ |3)^ ® \n p ,n t ) 

and I e^ v At the interaction time tint corresponding 
to the maxima of the gate fidelity in Fig. atoms are 
largely found in state \3) i ^\n p ,rit), and the relative 
populations and phase relations between the states of the 




FIG. 8: Average fidelity (top figure) and conditional phase 
shift (bottom figure) as a function of time for the three-level 
ladder scheme of Fig. and for N a = 10 8 , S p = IO7, 5 t = 
and g p = g t = O.OO227. The spontaneous emission rate is 
721 = 732 = 7 = 27r x 6 MHz. 



two photonic qubits are consistent with the "ideal" state 
of Eq. jnj. 



VI. EXPERIMENTAL VERIFICATION OF THE 
QPG OPERATION 

In this Section, two possible schemes for experimental 
implementation are discussed. First is the detection in 
the occupation number logical basis, and the second is 
the detection in the polarization logical basis. The two 
bases are identical from the point of view of theoretical 
treatment, however, their implementation is different in 
practice. 

Occupation number logical basis - In this section we 
describe a Michelson-like interferometer (see Fig for 
two-photons product state |1r)|1b) with the 'right' o~ 
circular polarization, where R (red) refers to the probe 
field and B (blue) to the trigger. We show that the inter- 
ferometer is able to reveal and measure the QPG phase 
shift. The probe and trigger fields with a bandwidth 
of 40 ~ 100 MHz (corresponding to 25-MO ns I/e half 
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width pulse duration |22() are separated in frequency by 
~7 THz and are resonant with the 87 Rb hyperfine tran- 
sitions DiF = 2 -»■ F' = 1 at 794.7 nm (377.228 THz), 
and D 2 F = 2 -> F' = 1 at 780.2 nm (384.225 THz), 
respectively. 



|1r>|1b) 



BS 



XPM 



FP! D B i 





FIG. 9: (Color online) Scheme of the proposed experiment 
for a measurement of the non-linear phase shift in a QPG. A 
Michelson-like interferometer with a two-photons input state 
|1r)|1b), probe and trigger, respectively, allows to measure 
the non-linear phase induced by the XPM on the logical ba- 
sis of the qubits, which coincide with the two lowest Fock 
states. Two intense classical fields at (tuner with Rabi fre- 
quency fii) and a c (coupler, with Rabi frequency Q4), are 
necessary to the 5-level XPM process. L is a lens for the 
mode matching in the EIT medium. BS a 50/50 beam split- 
ter. DL a delay line. FPi^ Fabry- Perot cavities. Dr.bi,b2 
are avalanche-photodiodes APDs. Left Inset - Scheme using 
a Sagnac interferometer for avoiding the optical path differ- 
ence. M is a mirror. Right Inset - Frequency spectrum [arb. 
units] for a 2 mm FP cavity length and finesse equal to 10 3 . 
The spectra of the probe and trigger photon are also shown 
in the plots as lines. 

The interferometer is realized with the help of a 50/50 
beam splitter (BS), using a Fabry-Perot cavities (FPx 2 ) 
instead of mirrors. The FPs reflect back the probe field, 
which is then superimposed on the BS and detected by an 
APD (Dr), and transmit the trigger field detected by an 
APD in each arm (Dbi e Db2)- This implies that only 
the trigger frequency is resonant with the FPs' cavity, 
which has a cavity length of 2 mm corresponding to a 
FSR of 74.85 GHz, while the probe frequency falls in the 
middle of the previous 93 th and 94 th FSRs. According 
to the photon bandwidth, a finesse of 10 3 determines a 
reflectivity for the probe field of 99.9% (see the spectra 
in the Right Inset of Fig. |5J. 

Since the frequency-bandwidths of the two photons are 
well distinguished and the FPs filter out the trigger field, 
this apparatus determines an interferometer for the probe 
field only. The coincidence probabilities P(R, Bl), be- 
tween Dr and Dbi, and -P(R, B2), between Dr and Db2, 
post-select the events in which the trigger photon is. In 



this case the coincidence probabilities are equal and given 
by 



P(R, Bl) = P(R, B2) = [1 + cos $] 



(21) 



where $ represents the phase difference due to the dif- 
ferent optical paths of the two arms experienced by the 
probe. In arm 2 a delay line (DL) is added to compen- 
sate the difference in the optical path and to scan for the 
interference pattern. In the Left Inset of Fig. 0a Sagnac- 
like version of the interferometer is shown, which allows 
for an auto-compensation of the optical path delay as the 
two arms coincide. 

When an EIT-based XPM system is considered in one 
arm, say arm 1, a non- linear contribution to the phase is 
added by the QPG, whether the trigger photon is present 
in the arm 1 or not. The XPM requires two intense 
classical fields resonant to the D 2 F = 1 — > F 1 = 1 
(384.232 THz) and D X F = 1 -> F' = 1 (377.235 THz), 
cr~ circularly polarized tuner field (at, fii), and a + cir- 
cularly polarized coupling field (a c , ^4), respectively. 
The coincidence probabilities P(R, Bl), between Dr and 
Dbi, and P(R, B2), between Dr and Db2, enable post- 
selection of the events in which the trigger photon is in 
the arm 2 and hence the interference pattern is given by 
Eq. (|21|1 . Instead a non-linear contribution to the phase 
is expected when detections on Dr are post selected by 
the detection on Db2- Fig. ^3 represents the diagram of 
the four amplitude probabilities after the action of the BS 
and the FPs on the two-photons state. The non-linear 
phase 4>n is added only to the first diagram. 
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FIG. 10: (Color online) Truth-table of the QPG for a log- 
ical basis of the qubits determined by the two lowest Fock 
states. On the right, diagrams corresponding to the probabil- 
ity amplitudes after the action of the BS and the FPs on the 
two-photons input state are shown. 

According to the truth-table and the amplitudes of 
Fig. the coincidence probabilities can be evaluated 
to be 

P(R,B1) = [1 + cos($ + 0)] /8 (22) 
P(R,B2) = [l + cos<f>]/8 (23) 

with $ = <f> + (0io — 20oo + 0+o) • 0+o is the phase due 
to the probe photon reflected in arm 1 with 'wrong' a + 
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circular polarization. The phases 0io and 0oo have been 
introduced in Eq. (|TJ as well as the conditional phase shift 
= 0ii — 0oi~0io+0oo of the QPG. The phase difference 
between the two interference patters determined by the 
two coincidence probabilities determines univocally the 
value of 0. In the case of an ideal QPG for which = it 
the two coincidence probabilities have opposite phases. 

Polarization logical basis - The previous proposal in- 
directly tests the QPG based on the XPM detecting the 
non-linear phase shift by a Michelson interferometer and 
coincidence measurements. A direct measurement of the 
truth table or a test on a general qubit state requires a 
control and measurement of a superposition of vacuum 
and one photon state. While the generation of a superpo- 
sition of vacuum and single photon state has been already 
achieved |22| . the measurement of such a su per position 
requires also a homodyne measurement |25l |26| . How- 
ever a logical basis for the qubits can be chosen as the 
circular polarization basis of the probe and trigger pho- 
tons. A test of the QPG will then require the detection 
of both photons thus avoiding a problematic production 
and detection of Fock state superpositions. 

The experimental setup is the same as in Fig. but 
the qubits are now encoded in the polarization of the in- 
put two-photon state. According to the truth-table and 
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FIG. 11: (Color online) Truth-table of the QPG for a logical 
basis of the qubits determined by orthogonal circular polar- 
ization basis. On the right the diagrams corresponding to the 
probability amplitudes after the action of the BS and the FPs 
on the two-photons input state. 

the amplitudes in Fig. ^] it is possible to derive the co- 
incidence probabilities P(Ri,Blj) and P(Ri,B2j), with 
i,j = {+,-}, as 



P(Ri,Blj) = [1 
P(Ri,B2j) = [1 



cos(O> + u .)]/8 (24) 
-cos$]/8 (25) 



with $ = <!>+ (0io — 20oo + 0(iei)o), where i®l is the sum 
mod 2. The phase is now given as cf)^ = 0ij — 0io — 0oj + 
0oo, where 0,o is the phase due to the EIT for the probe 
photon with polarization i and no trigger photon present, 
and same meaning for the other phases. Note that for 
i = j = ' = '1' one obtains the previous expression 
for the QPG phase in the logical Fock-state basis. Four 



possible choices of the probe and trigger polarizations 
determine the phases 



(26) 
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which satisfies the relation 

— |_ + ++ = (f>. In the case of an ideal EIT 

!> +0 = 00+ = 0++ = 000, 0-+ = 0-0 = 4>R 

- 0o- = 0s, where 4>a, b are the phases ac- 
quired by the single photons, we have = 0, and 

|_ = <p^ = ++ = 0. In a way the phases between 

the two coincidence interference patterns allow a mea- 
surement of the QPG phases in the diagonal basis, i.e. in 
the single qubit states for which the only non zero phase 
shift is the conditional phase shift of the QPG. 
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FIG. 12: (Color online) Scheme of the proposed experiment 
for a complete characterization of the QPG. Two photons in 
the a + polarization state (logical state \0r) <%> \0b)) are trans- 
formed by a QWP, corresponding to two Hadamard single- 
qubit gates, and then to the XPM medium (QPG). A Fabry- 
Perot cavity with the same parameter as before transmits the 
trigger photon and reflect the probe to two tomographic mea- 
surement systems (Tr : b) and detected by APDs. 

General polarization qubit input states - The polariza- 
tion logical basis allows for a direct observation of coher- 
ence and production of entanglement as necessary condi- 
tions for a QPG. Assume that the information is encoded 
in the polarization state of two photons, and then sent 
into the EIT-based XPM system for QPG, as shown in 
Fig. El The output photons, as shown in Fig. are 
split by a dichroic mirror (a tilted FP cavity with the 
same parameters as before) and then collected in two 
APDs (D^ and T>b) for coincidence counting. In front of 
each detector a tomographic system |27| constituted by 
a QWP, HWP and a PBS, is placed for the complete re- 
construction of the polarization state of the output, thus 
providing the information on coherence properties of the 
gate. It has also been shown Q that an input state for 
the QPG given by [(|a+) + |a^)) ® (|<r+) + |<t b )))]/2 can 
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quantify the entanglement of the output state, for which 
the CHSH inequality is 2\J\ + sin 2 <fi where 2 is the upper 
classical limit. 

VII. CONCLUSION 

In conclusion, our study shows that the implementa- 
tion of efficient EIT-based nonlinear two-qubit gates for 
travelling single-photons is possible, even though experi- 
mentally challenging. The main limitation is due to the 
existence of a trade-off between the size of the CPS and 
the fidelity of the gate, limiting the achievable gate fi- 
delity in the stationary regime, but which can be par- 
tially bypassed in the transient regime. Since this trade- 
off is a general consequence of the coherent interaction 
between the atomic medium and the single-photon wave- 
packets, we expect that these considerations apply to all 
EIT-based crossed-Kerr schemes 0, , regardless the 



specific level scheme considered. Instead, this consider- 
ation does not apply to situations where the nonlincar- 
ity comes from independent processes (e.g. collisions of 
dipole-dipole interactions) [14( , nor the similar solid-state 
based processes [24). 
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